The evaluation of directionally writhing polymers. 
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Abstract. We discuss the appropriate techniques for modelling the geometry 
of open ended elastic polymer molecules. The molecule is assumed to have fixed 
endpoints on a boundary surface. In particular we discuss the concept of the 
winding number, a directional measure of the linking of two curves, which can 
be shown to be invariant to the set of continuous deformations vanishing at the 
polymer's end-point and which forbid it from passing through itself. This measure 
is shown to be the appropriate constraint required to evaluate the geometrical 
properties of a constrained DNA molecule. Using the net winding measure we 
define a model of an open ended constrained DNA molecule which combines 
the necessary constraint of self-avoidance with being analytically tractable. This 
model builds upon the local models of Bouchiat and Mczard (2000). In particular, 
we present a new derivation of the polar writhe expression, which detects both the 
local winding of the curve and non local winding between different sections of the 
curve. We then show that this expression correctly tracks the net twisting of a 
DNA molecule subject to rotation at the endpoints, unlike other definitions used 
in the literature. 



PACS numbers: 02.40.-k, 82.35.Lr, 87. 15. ad, 02.40.-k 
1. Introduction 

The linking of a pair of space curves has been quantifiable since Gauss, who 
developed an expression for evaluating the integer degree to which the paths of 
two celestial bodies were linked (Epple [13] )• In contemporary study this measure 
has been widely used to characterize the helical inter-linking of the two phosphate 
backbones of the DNA molecule (see Bates and Maxwell [2] for an overview of the 
role played by topology in the analysis of DNA) . The geometrical study of DNA has 
been greatly assisted by a theorem of Calugareanu (Calugareanu [10], [IT] . White 




Figure 1. An example of supercoiling, of a section of DNA. The two figures 
are representations of a section of DNA molecule. In (a) the DNA's axis coils 
around itself to form a loop type structure. This loop is said to have writhing. 
This writhing is not present in (b). The two figures are interchangeable by an 
appropriate set of deformations applied to the axis. This figure is reprinted from 
Travers and Muskhelishvili [32 . 



[35] Fuller 1971 [15]) which separates the linking into components, representing the 
self linking of the helix's central axis (the writhe W) and the number of times 
the backbones wind about this axis (the twist T). Denoting the linking as C 
Calugareanu 's theorem gives us 

C = T + W. (1) 

We must state at this point that this theorem was defined for closed pairings (e.g., 
DNA Plasmids) and is not, in general, well defined when the constituent curves are 
open-ended (though in certain cases it can be shown to hold, van der Heijden et al 
|33j). Writhing has a special place in the study of DNA as it is used to quantify 
the phenomenon of supercoiling, in which the DNA molecules axis entangles with 
itself such that long DNA molecules can occupy a small area in the structure of a 
cell (Figure [TJ. 

If we consider a space curve x(£) (where t is and arbitrary parametrization) , 
which could represent the helical axis of a DNA molecule, its writhe W can be 
quantified as 

*II |x(i) - x(OI 3 

where T x (s) is its the unit tangent vector of x. Note that this expression is 
independent of parametrization. One can interpret Q as a directional average 
as follows. Consider a Cartesian direction b representing a viewing direction. 
We project the curve onto a plane perpendicular to this viewpoint as shown in 
Figure |2](a), with overlapping sections shown drawn above or below each other as 
demonstrated in the figure. If the curve is oriented we can apply a sign to each 
crossing as depicted in Figure [2j The sum of these signs defines a measure of the 
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Figure 2. The diagram (a) represents the planar projection of a single space 
curve along a direction b. Diagram (b) is the result of this projection as viewed 
along the direction of projection. Diagrams (c) and (d) represent example planar 
writhing calculations, using the crossing rules shown to the right. 




Figure 3. A looped section of curve is split into two sections by its maximal 
point along z. These sections are labelled Xi and X2. The two sections of curve 
which are separated by the local maximum in z, can be seen to wind around each 
other in the x-y plane. Such contributions to the writhing geometry of the curve 
are non local and are ignored by the Fuller writhing expressions. 

planar writhing for a particular projection (see examples in Figure [2]^c, d)). The W, 
as defined by g can be shown to be the average over all directions of projection 
(see Pohl 1968 [23J or Aldinger et al jl]). As a result of this non-locality its use 
in analytic calculations can be limited (Bereton and Shah [7]). In 1978 Fuller [16] 
introduced two single integral expressions for evaluating the writhing of closed space 
curves, which are defined in terms of the rotation of an orthonormal trihedron whose 
orientation is defined at all points by the local geometry of the curve x. However, 
both expressions can only accurately evaluate the writhing of the axis mod 2. This 




Figure 4. A depiction of optical tweezer experiments. The molecules can 
be manipulated by optical tweezers. The magnets are used to manipulate the 
paramagnetic bead. The tow applied forces are a Torque force imparted by 
rotating the bead and s vertical force stretching the molecule. The right figure is 
a plectonemic structure which results form the application of a significant Torque. 
This figure is reprinted from Bustamante et al [8]. 

is due to their ignorance of non-local (global) writhing (see figure [3]). That said, 
if certain necessary, but rather geometrically complex, conditions are met they can 
be used for both analytic calculation and efficient numerical evaluation (Starostin 
[30] . van der Heijden et al [33]). We do not discuss these conditions here as the 
aim is to provide a framework for evaluating the geometry of curve for which such 
restrictions are not necessary. 

The development of magnetic tweezer and optical trap techniques have allowed 
researchers to investigate the in vitro response of DNA molecules to torsional stress 
(see Smith et al [31], Bustamante et al [8] and Bustamante et al for a review [9]). 
Of particular interest has been experiments in which an open ended section of a 
DNA molecule is held fixed at one end and manipulated at its second end by a 
paramagnetic bead (see Figure [4]). The molecule is subject to both a stretching 
force along this vertical direction and a torque perpendicular to the stretch force. 
We emphasize at this point that the forces controlling the molecule are either along 
or perpendicular to a specific Cartesian direction. Attempts to match the findings of 
these experiments have been the concern of a number of theoretical papers. These 
models focus on a polymer chain and require a partition function evaluation over 
the set of allowed geometrical configurations available to the molecule. In the limit 
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Figure 5. A depiction of typical set-up used to manipulate a DNA molecule 
in experiments. The molecule depicted as two phosphate backbones wound into 
a helical shape, typical of DNA molecules, is attached to a fixed surface and a 
paramagnetic bead. The bead can be rotated (about a specific axis) by a magnetic 
field in order to apply a torque to the DNA molecule. 



of high force the molecules elastic behaviour is well understood using the worm- 
like chain model (Bustamante et al [8], Vologodskii |34j). Here the molecular axis 
only deviates slightly from pointing along z (the deviations are a result of random 
thermal fluctuations). There is currently more interest in the so called entropic 
elastic domain in which the set of allowed configurations can involve a notable 
degree of writhing (Vologodskii and Marko [TH], Bouchiat and Mezard [5j, Moroz 
and Nelson [20], Sinha [26]). In particular Bouchiat and Mezard comprehensively 
detailed a model [|J which was shown to provide a good fit to the experimental data 
in this domain (similar results were also obtained by Moroz and Nelson |20j). 

The aim in such statistical models is to identify all possible energy partitions (a 
continuous set in such macroscopic models, see Bouchiat and Mezard [5]). In order 
to do this Bouchiat and Mezard devised an expression for the axial writhing of the 
molecule in terms of Euler angle rotations. One interpretation of this measure is 

| Termed the Rod like chain (RLC) model by the authors, but also referred to as the south avoiding 
worm like chain in Samuel et al |27j . 
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an application of the local Fuller writhing expressions, though in a different form 
from that originally specified by Fuller Q16J). It defines the local writhing of the 
molecule's axis about a specific direction (in their paper it is the vertical direction 
z). In Cartesian form the expression (which we here denote W z ) is the rate at which 
an orthogonal frame rotates about z: 



A requirement, resulting from the ill definition of the Euler angles at 6 = it, is 
that the axis tangent curve is restricted from pointing along —i anywhere along 
is length. It appears to be the authors' claim that this restriction allows the local 
•writhing expression to replicate the behaviour of the better understood non-local W 
expression (J2|. In doing so the claim is made that this model of south-avoidance, 
that is the molecule's axis cannot point vertically downwards at any point along 
its length, will accurately capture the geometry of self avoidance (the appropriate 
physical obstruction to the systems global geometry). It is this particular issue 
which has led to a significant degree of debate as to the model's validity. Several 
papers, framed from a geometrical view point, have demonstrated that the local 
writhing expression does not accurately replicate the behaviour of the total writhing 
expression ^ (Rossetto and Maggs [21] and in particular Neukirch and Starostin 
[2Tj and [22])- On the other hand Samuel et al [21] and [2E] take the view that 
the configurations, for which (|3| is in error, offer a negligible contribution to the 
full partition space, thus explaining the apparent success of Bouchiat and Mezard's 
model. We intend to show that this discussion is unnecessary as a more suitable 
writhe expression can be employed for open-ended, restricted polymers. 

This note contains the derivation the path integral based model for defining and 
evaluating the appropriate non-local partition function of an open ended polymer, 
such as the DNA molecules studied in the micro-manipulation experiments we have 
discussed. With regards to the experiments, a key factor in the modelling of the 
molecule is the directional nature of the applied forces. The forces applied to the 
field are transmitted through control of the paramagnetic bead (figure [5]) to which 
the upper end of the molecule is attached. This bead is placed in between the north 
and south poles of a magnet. By altering the polarity of this magnet the bead 
can be rotated by applying a torque force to the molecule. This force is always 
perpendicular to a chosen direction. In this note we shall chose this direction to 
be vertical (z) following the convention used by Fain et al [14] and Bouchiat and 
Mezard [5]. That is to say the bead rotates in the x-y plane. Secondly a force F 
can be applied to the bead stretching it along z. 




(3) 
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The paper will take the following structure. In section [2] we discuss the energy 
expression of an elastic ribbon representing the polymer molecule. The expression 
is derived by considering the rates of rotation of an orthonormal frame attached 
to a curve representing the polymer's axis. In particular the twisting expression 
differs from that in the literature ([H],[S],[2n])- However, it is shown that despite 
this difference the energy expression for the ribbon is the same as that used in the 
cited literature. 

In section [3] we introduce the net winding expression. This is a measure of the 
directional linking of a pair of curves. A discussion relating this expression to the 
confined polymer experiments we have discussed in the introduction demonstrates 
that the net winding is the correct geometrical constraint, which should be used to 
restrict the allowed configurations of the ribbon. 

In section [4] we detail the correct writhing expression for confined polymers. 
This measure was previously introduced using the net winding measure and is 
termed the polar writhe (Berger and Prior |4J), however, it is derived here in a 
different manner. We use the correct form of a magnetic monopole field in order 
to derive the polar writhe. This derivation supplies an interpretation for the link 
between the quantum mechanical problem of a symmetric top in a magnetic field 
and a constrained polymer molecule, subject to torsional and stretching forces. This 
link was first suggested by Bouchiat and Mezard without a geometrical justification 
which we present here. Further to this we use two example curve studies, a helix 
and a deformed spacecurve whose end points are fixed, in order to demonstrate that 
the polar writhe measure has the required geometrical properties for constrained 
polymers. 

In section [5] we derive a partition function, using the net winding measure, 
which satisfies the non local self-crossing restriction, along with the possibility 
of being evaluated analytically. This is achieved by splitting the net winding 
measure into local and non local components. The local contribution is decomposed 
into directional twisting and winding contributions in a manner similar to that of 
Bouchiat and Mezard (though with differing twist and writhe expressions). As the 
energy expression is defined in terms of local geometry we use the non local winding 
contribution as a constraint which enforces the non locality of the model. 

2. Defining the energy of the polymer 

2.1. Geometrical background 

In what follows a space curve x(t) shall represent a three-dimensional vector 
(xj., x y , x,), depending continuously on an arbitrary parameter t, for t G [a, b}. The 
Euler angles in this note will be the set (0 , 6, ip) with a rotation in the x-y plane, 
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6 a rotation in the y-z plane for z and ip a rotation in the x-y plane, following the 
rotation 9. In addition to the Cartesian system and the Euler angle set, defined in 



section Appendix A, we will use the spherical polar coordinate system (r, 9, </>), with 
theta the polar angle and <p the azimuthal angle. It shall be made clear when we are 
discussing spherical polar or Euler angle parametrisations. All such curves will be 
considered smooth, where smoothness implies x is at least C 3 differentiable for all t. 
Differentials are denoted with a prime (ire x'(t) = dx/dt). 

2.2. Ribbons 

The ribbon is a mathematical construction used to represent polymers such a DNA 
molecules (Fuller [IS])- Consider a space curve x(t) and a second curve y(i), also 
parametrized by t, such that y(t) = x(i) + eu(t), where u(t) is a vector normal to 
T(t) for t £ [a, b\. This will naturally wrap itself around x as shown in Figure [6 If 
e is sufficiently small (usually 1) we can assume that y is disjoint from x, that 
is x and y never cross (Hirsch [IB]). Such a construction will be denoted i?(x, y). 
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2.3. The enegry components 

In order to create our partition function we must define the energy required for the 
polymer to assume each possible configuration. In Bouchiat and Mezard the elastic 
energy of the polymer molecule was define using a set of rotations (parametrized 
using the Euler angles) of an orthonormal frame ^T, U, . The rotations were 
relative to a fixed rectilinear reference frame pointing vertically (along z). Let 
T z = T • z be the vertical component of the tangent vector, and 9 = cos" 1 T z . Also 
let a give the sign of T z , i.e. 



1, if < 6? < tt/2, 

-1 if tt/2 < 6 < 7i. U 



We will multiply our rotations by a. A physical justification is given in Appendix 
A. We shall see that, whilst this does affect the expression for the net twisting of 
the molecule, it does not affect the energy expression using in Fain et al [13] or 
Bouchiat and Mezard [5]. However, we shall see later in section [5] that the difference 
in twisting does have an effect when we apply the necessary geometrical bounds on 
the set of configurations accessible to the molecule. 

We define the angular velocity (ft) of our orthonormal frame as 

fl = uiT + w 2 U + w 3 V, (5) 

where the u values are scalars representing the rates of rotation about that particular 
direction. For the sake of brevity we do not detail the derivations. Bouchiat and 
Mezard have already detailed the procedure for deriving the three angular velocity 
rates. We have repeated this procedure whilst applying the rule set defined in Q 
and find the following expressions 

0J\ = h costal — . (6) 

az az 

U2 = Vi — . (7) 

az 

lu 3 = Oi sin 6^-. (8) 

az 

2.4- The energy expression 

In general the energy of a particular configuration of an elastic polymer, modelled 
as a ribbon, has three components. The first component results from the bending. 
Following the work of Fain et al [Hj and Bouchiat and Mezard [5] this can be 
quantified by the two angular velocity components U2 and u 3 which represent the 



The evaluation of directionally writhing polymers. 10 
rate of rotation of the tangent vector about U and V, leading to the following 
density, 



2 



where A is the bending rigidity constant. The second contribution comes from the 
twisting of the polymer's edges about its central axis. This can be defined by the 
density 

C 2 C ( dtfj nl d<f)\ 2 
E' twlst = -ujl = -[o l -^ + \ C ose\^) , (10) 



2 1 2 V dz 1 1 dz / 

C (di) „d<A 2 
= — -- + cos6» -- . (11) 
2\dz dz) K J 

Here C is the twist rigidity constant. Finally to complete our expression we allow 

for the molecule to be stretched. In this experiment the molecule is stretched along 

the z axis. A uniform stretching force (F = F ■ z) is applied at the free end of the 

molecule. This is a potential and acts at all points along the molecule based on 

their orientation. If the molecule is moving upwards in z then the force will have a 

negative value. If it is moving downwards it will have a positive value, giving 

Fcosfl 

^stretch ~ ^ rp ■ \ lz 7 

3. An invariant linking expression — The net winding number 

Having discussed the energy expression for our cylindrically symmetric elastic 
polymer we must now consider the geometrical restrictions placed on the molecule. 
It is at this juncture that our work begins to differ significantly from other treatments 
of the problem (Fain et al [H], Bouchiat and Mezard [5], Moroz and Nelson [20J, 
Rossetto and Maggs [23] > Samuel et al |2"7]). For this discussion we shall assume 
the torque and stretching forces take on fixed values (though we assume they are 
arbitrarily chosen, such that this discussion is generalized). In this scenario we 
would like a measure of the polymers geometry, which remains fixed for fixed torque 
and stretching values. In Berger and Prior [1] just such an invariant was defined, 
which we refer to as the net winding C. It represents the extent to which the axis of 
the polymer and a second curve, (possibly representing one of the protein strands 
of a DNA molecule), wind about each other perpendicular to z (see figure [7]). It is 
of note that this expression is not in general equal to the Gauss linking number (£), 
which is a well known topological invariant for closed polymer ribbons. Authors 
have suggested artificial extensions of the polymer-ribbon which allow the use of the 
Gauss linking invariant (Rossetto and Maggs [25], Samuel et al [27]), we show in 
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Figure 7. Two sections of curves x and y occupying the same z range. Their 
linking, along z, can be defined in terms of the vector r(z) joining two sections of 
curve x(z) and y(z) in the x — y plane. Also depicted is the angle 9(z) which r 
makes with the x-axis. is a combination of the two Euler angles <f> and ip. 



this section that this extension is not necessary, as the appropriate invariant can be 
described without the need for an artificial extension. This measure will be split 
into two distinct contributions, the local winding Ci and the non-local winding C n i. 

3.1. The local winding 

Consider a section of the polymer ribbon i?(x, y), parametrized by z, for which 
both sections are travelling upwards in in z (x.'(z) > and y'(z) > 0), over a range 
z £ [zmin, z max]- We define a vector r, originating at x which will lie in the x-y plane 
such that its tip lies on y(z) for all z G [z min , z max ], 



(see Figure [7]). We wish to measure the rate of rotation of this vector in the x-y 
plane, which can be used to represent the extent to which the two curves wind about 
each other, perpendicular to z, over this range. In order to do so we must evaluate 
the rate of rotation of the vector as we increase z, in Euler angle rotations this would 



By integrating over z we obtain the net winding of the ribbon R(x, y) perpendicular 
to z. This result was obtained by Bouchiat and Mezard in a different fashion. 
The same expression applies as for sections of the polymer which are travelling 



r( z ) = y(z) - x(z) 



(13) 



be 




(14) 




Figure 8. A depiction of a section of ribbon i?(x, y) which has both local winding 
and non local winding. Both curves are split into three sections (x 1; x 2 , x 3 , y 1; 
y2 and by their turning points along the z-axis. Shown on the diagram are 
two planes z = z m i n and z = z max which bound the range of z values shared by 
X! and y 2 . Also depicted are the directions of r 12 at z = z min and z = z max . we 
used these vectors to evaluate the non-local winding of the ribbon. 



downwards in z (x'(z) > and y'(z) > 0). We have termed this expression the local 
winding Ci as all contributions arise form the winding of curve sections Xj and y~j 
belonging to the same section of the ribbon. 

3.2. The non-local winding 

In order to render this measure a topological invariant (to the set of deformations 



vanishing at its endpoints, see section 3.3) we must consider a further set of 



contributions due to the winding of the polymer between distinct sections of the 
polymer. Consider a polymer configuration which has two turning points in the 
z direction (Figure [8]). There will be additional linkage of the two curves for the 
range of z values which share a mutual z range as marked on Figure [8} Consider for 
example the sections marked xi and y~2. In order to define the net winding of this 
pair we use a vector r 12 (see Figure [8l which is defined as 



ri2(z) = y 2 (z) -^(z). (15) 
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By defining the infinitesimal rate of rotation as described above we can recover the 
net winding between these two sections by integrating over their mutual z range 

rain y raa 
12 '^12 



•y d \ fTnm ^raax 
Z t yz 12 , -6 



£ 12 (R, ,5") = I I'"" - + ) di, (16) 



' -/nu n \ dz dz 

" z 12 v 

where the labelling and <pi2 indicate that they represent the required Euler 
rotations necessary for the top-vector combination to trace out the polymer section 
i2(xi,y 2 ). The minus sign occurs as the two sections are travelling upwards and 
downwards along z respectively (see Figure [sj) . We can generalize this construction 
to cover all geometries which are possible for a polymer bound at two specific points, 
which are distinct (this is exactly the case for the DNA molecule discussed above). 
For a polymer constructed from curves x and y which have n — 1 and m — 1 turning 
points along z we consider the set of vectors attached to the top centred on 
section x« which points to section y^. The total net winding of a polymer R(pc,y) 
, about z between two planes z min and z max , is then the sum of all contributions 
between sections Xj and yj, 

n m 

/2(R, Z m i n , Z max ) ^ ^-"iji (1^0 



EE 



max , 



+ ^)dz 



V dz dz 

i=l j=l *ij x 

Here the local contributions will be those for which i = j, which we have labelled 
Ci. We call the contributions for which i ^ j the non-local contributions, we shall 
label these L n \. A key factor in our final model will be the distinction between the 
local and non-local contributions. 



3.3. Topological invariance and directionality 

We can show that the net winding has a restricted sense of topological invariance, 
which is sufficient for open-ended ribbons whose endpoints are constrained. Before 
we discuss this theorem we introduce the concept of the base state. Fain et al [14] 
show that the minimum energy state for such constrained polymers are braided 
vertical lines such as that marked (b) Figure [9] We shall call such states the base 
state. We assume that the curve can be arbitrarily deformed within some region of 
space, where the endpoints attach to the boundary of the region. What we shall 
show in the following theorem is that writhed configurations such as (c) (Figure [9]) 
which result from deforming the base state (b), subject to its end orientations being 
fixed, will have the same net winding. 




Figure 9. Depicted are three possible configurations of the constrained polymer 
(representing a DNA molecule for example). The configuration depicted in (a) is 
a rectilinear configuration which would only be present in the case in which there 
have been no turns of the bead and not linking of the polymer strands before 
attachment. The configuration depicted in (b) is a braided configuration which 
would result from twisting the bead as shown in (a). The number of windings 
of the two strands of the molecule about each other would be proportional to 
the number of turns applied to the bead (assuming we started with a rectilinear 
configuration). The molecule's axis is a straight line for (b) so it has no writhing. 
The figure in (c) is a configuration in which the number of turns applied to the 
bead is the same as (b). It has both non-local and local windings. In (c) the axis 
itself has a non-trivial geometry. Both configurations can be shown to have the 
same net winding. 



In Berger and Prior the following theorem was proven. We restate it here such 
that we can use it in the following discussion to justify the use of the net winding 
as the correct topological restriction for constrained polymer molecules. 

Theorem 1 Consider an continuous open ended ribbon i?(x, y) (both x and 
y are assumed to be continuous). We shall assume this ribbon is of at least C 3 
differentiability. Let z m \ n be the height of the lower endpoint and z max the height of 
the upper endpoint. We restrict the motion of the ribbon to z min < z < z max . We 
define the set of restricted ambient isotopies, which vanish outside z m i n < z < z ma , x , 
as the end-restricted ambient isotopies. Then the net winding, defined by 
equation (18) is invariant to all end-restricted ambient isotopies. 

For proof see Berger and Prior j3], the pertinent theorem in that paper is 
theorem 2, which in turn relies on theorem 1 (theorem 1 of [3J, not the above) . 
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Briefly the net winding of a ribbon of infinite length, is shown to be equal to the 
Gauss linking number (theorem 1 of [4]). Consider a section of this curve (the section 
of interest) to which the deformations are constrained, and at whose end points the 
motions vanish. The net winding measure can be shown to be equal to the linking 
of the infinite assembly minus the net winding of the remaining ribbon sections (i.e. 
everything but the region of interest) which must be invariant as all motions vanish 
outside the region of interest. Thus the net winding of the section of interest is also 
invariant. 

It can be seen that, in terms of the systems physics, the net winding is the 
appropriate measure for evaluating the linking of the constrained DNA molecule. 
Consider a rectilinear polymer configuration as depicted in Figure |9^a). We then 
apply a Torque in the x-y plane as in the experiments (again see Figure [9|a)). 
This torque will cause the polymer to form the base-state braided structures (see 
figure |9](b)) with the number of turns of the braided structure corresponding to 
the the number of turns applied to the bead, this would naturally be measured 
by the net winding (18). It is possible, given a significant number of turns of the 



bead, that the polymer structure can buckle to form structures such as Figure |9[c). 
These structures have significant writhing (we shall detail the appropriate expression 
for evaluating this writhing in section [4]). What should be noted is that we have 
demonstrated in theorem 1 that the net winding will stay fixed, assuming the end 
points are prevented form rotating; the winding is merely distributed in a different 
manner. It is for this reason that we use the net winding as means of con straining 
geometrical configurations available to the system, for a given number of turns of the 
paramagnetic bead. This suitability derives entirely form the directionally specific 
nature of both of the forces applied in the experiments. We must reiterate that the 
net winding is not, in general, equal to better known Gauss linking expression L 
which is defined as an average over all viewpoints and it is physically appropriate, in 
this case, to use the net winding rather than the C This fact was already recognized 
by Bouchiat and Mezard who, derived their own directional winding expression. We 
have simply defined the appropriate form of the required directional invariant. 



4. Detecting the writhe using magnetic fields 

In this section we define the appropriate form of directional writhing of an open 
ended polymer's axis. In section [5] this expression will be used in order to define an 
analytically implementable partition function which incorporates the self avoidance 
condition. This writhing expression, termed the polar writhe, has already been 
specified by Berger and Prior |4], using the net winding expression as a basis 
for construction. The paper was primarily written from a purely geometrical 
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perspective. We choose to introduce this expression here in a different fashion, 
using the language of magnetic fields. Bouchiat and Mezard derived a version 
of the directional writhing expression which was shown to be equivalent to a line 
integral over the vector potential of a field generated by a thin magnetic solenoid 
which extends from infinity to some endpoint x (in effect this is the magnetic 
field of a Dirac monopole placed at position x). They used this fact to develop a 
Hamiltonian expression based on a quantum symmetric top subjected to a magnetic 
field, though there was no real suggested physical interpretation for this link. In what 
follows we shall define the appropriate form of the directional writhing, utilizing 
the correspondence between electromagnetism and the geometrical language of fiber 
bundles first demonstrated by Wu and Yang [36] . In doing so we will assign a physical 
interpretation to this equivalence, with regards to the DNA-polymer experiments 
we are discussing. 

We shall construct our writhing expression using a tethered magnetic monopole 
(of charge q = 1) and its vector potential. As magnetic monopoles do not possess 
globally defined vector potentials (and, as far as we know, do not exist in nature) 
we imagine a thin bundle of magnetic flux $ = q/(4ir) extending from z = — oo to 
the point x, confined to a thin vertical tube. At x the flux escapes from the tube 
and spreads equally in all directions, mimicking a magnetic charge. Our set up will 
be equivalent to the rotated symmetric top description described in Appendix A 
(section Appendix A), with the charge replacing the top and the vector potential 
the attached directing vector. We consider a curve x which represents the axis of 
our polymer. We shall first consider a section of polymer moving only upwards along 
z. We start the magnetic monopole at the base of the polymers axis (minimum z 
value) and attach the following vector potential 

|a| (|a| + z • a) ^ 
Here a represents a vector (x* — x), with x being a point belonging to the field line 
and x' any other point between the planes. This vector potential represents one 
possible vector potential used to describe the Dirac magnetic monopole (Dirac 1931 
[12J). This potential is well defined for all x* except points vertically below x. For 
our purpose we shall demand that a is of unit length, further, we shall demand that 
for each x the point x* will be such that a points along the direction of the unit 
tangent T giving 

A N = " XT ^ (20) 
1 + z-T 

We now consider the potential at a single point x and consider an infinitesimal 
change in the geometry of the field line dT. This will mean the field defined by the 
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magnetic monopole will deform. We take the dot product of the potential and this 
infinitesimal change to define the rate of rotation of this potential 

■ dT = • dT. (21) 

(l + i-T) 

A key geometrical property of the potential A N is; whichever direction the field 
line is pointing the potential will always be projected into the x-y plane. Thus the 
winding is defined in the x-y plane, as we would require. As the field line x always 
has a positive gradient along z we can state that for each z value the potential 
vector will be unique, that is none of the potential vectors will intersect. Using 



dT = T(z)'cLz we can rewrite (21) as an infinitesimal in terms of z 



A n . dT = ^( TxT V (22) 

(i + i-T) 

In spherical coordinates this gives 

A N -dT= (l-cos0)^dz. (23) 

This expression is the same in spherical coordinates as its would be in Euler angles. 
The the third rotation if> would be assigned to rate at which the vector potential 



rotates about the axis (see section 4.2) 



We can interpret this geometrically. Consider an observer in the moving 
reference frame of the monopole as it moves up the field line (up in terms of z). 
This observer would note a continuous change in the shape of the potential field 
generated by the flux point source (assuming x is not planar). The change in the 
potential field is due entirely to the geometry of the polymer's axis. It is the change 
in this potential which shall be used to evaluate the writhing of this curve section. A 
means by which this observer could track the the fields changing geometry would be 
to track a point whose relationship to the observer is fixed, in this case the tangent 



vector. The expression (23) represents the rate of continuous change of the field as 
the magnetic monopole is transported along the axis x. We integrate over the range 
z £ [zmin, z max] in order to obtain the total change in the axial geometry, we divide 
this by 27r to obtain a measure of the writhing of the section x 

1 f Zmax d(f) 
W p/ (x, z min , z max ) = — (1 - cos6>) — dz. (24) 

2n L . dz 



This is exactly the directional writhing expression used by both Fain et al 
and Bouchiat and M ezard [5] to model the axial geometry of the constrained DNA 
molecule. In particular Bouchiat and Mezard note the issue of the singularity and 
define their partition function model such that the direction — z is avoided using a 
cut off potential. 
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We take a different route to handling the singularity, we choose to define an 
alternative potential for sections of curve with a negative z gradient. We consider 
a section of curve x, again representing the polymers axis and whose gradient is 
negative along z. In this case we begin the monopole at the maximum z value in 
the range in which this curve section is defined. The vector potential we attach to 
our monopole will be 

. c Z X cl , . 

|a| (|a| — z ■ aj 

Again set a = T and take the dot product with an infinitesimal change in the 
tangent direction 

X s ■ dT = • dT. (26) 

(1-5- T) 

In spherical coordinates we have 

A s ■ dT = - (1 + cos 9) ^dz. (27) 

dz 

We integrate this over the range z G [z max , z min ] in order to recover the directional 
writhing of this section. 

Wpi(x, z min , z max ) = — / A 5 -dT = — / (1 + cosfl) -^dz.(2S) 



Where we have reversed the integral's limits on the far right of (28) 



4-1. The non local contribution to the writhing 

It was shown by Wu and Yang 1975 [36] that we can avoid the singular 
representations by ascribing the potential A N for sections of x with positive 
gradients along z and A s for sections whose gradient along z are negative. In 
order to make this system consistent they show it is necessary to apply a gauge 
transformation in regions of overlap. In this case we must consider regions which 
share mutual z ranges, we must define a gauge transformation for each z value, 
that is 9 = tt/2. For such regions we employ the gauge transformation (in spherical 
coordinates) , 

A N - A 5 = WA eq = 2^dz, (29) 

az 

What meaning could be ascribed to a line integral the gauge transformation VA eq 
for a range of overlap z G [z m im z m ax}7 Consider a general curve which has sections 
moving both positively and negatively about z. This will naturally mean there 
are sections of curve which wind about each other about the z axis. Consider two 
sections of x, Xj and x 3 - which share this mutual range z G [zi,z 2 ]. We define a 
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vector Yij{z) = x 3 -(z) — Xj(z) and denote its orientation in the x-y plane as Q(z) 
(the orientation from the x axis). This is of course the vector used to define the 



net winding (18), excepting that the two curve sections are part of the same curve 
x, indeed we will be defining the self winding of these two sections. The rate of 
rotation of r is 



VA- ■ r = Z - Ti * Z) f X **( Z) dz (30) 
viz) r 



thus, in terms of Euler angles, we have 



22 Z ■ Yij{z) X Tij\z) __ f Z2 ( dlpij d</>„ 



r„M ? - Jz, V d, + -Jf) A ~- (31) 
As with the net winding we consider the contribution due to ijj which will be 
identical. This integral gives the net angle through which the vector potential z x re- 
winds, about the z axis. We call this term the non local polar writhing W pn i as is 
defined in terms of a vector potential linking distinct sections of the curve x. 

4-2. The Polar writhe 

Consider the axis of a ribbon, represented by the curve x, which spans a range 
z £ [z m in, z max] and which has n — 1 turning points along z (4^ = 0), and which can 
be split into n sections. We can define the net local polar writhing as 

W p z(x, z min ,z max ) = — / (1 - \cos9i\)-r^dz. (32) 

i 

We further define the net non-local polar writhing, W pn i. Consider sections of the 
curve Xj and Xj, which share a mutual range of z values z G [z^ m , z^ ax ], we have 

T x ( dipij (\o,j 



w ,(x* ■ z ) = rv^ f" ( 

" ypnl ^min i <~"max ) / 7 o / \ 

8=1 j=l J Z ij v 



dz dz 



+ 7 dz. (33) 



Finally we define the polar writhing W p , of the curve, as the sum of these two 
contributions 

^^p(x, z m i n: z max ) Wpi(x, z m i n: z max } -\- Wp n ^(x, z m i n , z max ^j. (34) 

4-3. A comparison the the writhe result of Bouchiat and Mezard and Fain et al 

We can compare this result to the expression used by Bouchiat and Mezard in order 
to understand the various components of the polar writhe measure. Consider first a 
configuration whose gradient along z is always positive. In this case we would only 
need to apply the vector potential A N and the polar writhe is defined entirely using 
the vector potential A N . In effect the tip of A N would draw out an imaginary curve 
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y creating a ribbon i?(x, y). Such a ribbon would have C and twisting measures. 
We can define the local twisting as the integral over the rate oj\ (|6|) divided by 2n 



1 [ Z ™ax / fajj d(j)\ 

T(R, z min , z max ) = — / I — + | cos 9\ — j dz. (35) 

We note the difference between the net winding (here given entirely by its local 
contributions) and the twisting gives the required polar writhing expression in this 
case, that is 

^i^Ttj z m i ni z max ) T~( y R ) z m i ni z max ) Wp/(.R, z m i n , z max ) . (36) 

This result was also obtained by Bouchiat and Mezard. However, differences arise 
for polymer configurations for which there are downward travelling components. 
We note that, due to the | cos^l in our twisting expression (36) is also true of the 
A s contributions to W p ;. This is the first major difference between this directional 
writhing formulation in contrast with the work of Bouchiat and Mezard. In this 
case northern and southern local writhing are given equal weighting. 

The second difference results from the contribution to the writhing resulting 
from the gauge transformation VA e? . A polymer with at least one turning point 
along z will have a non local writhing contribution. Rossetto and Maggs [25J 
and Neukirch and Starostin j2T] have already noted that curves with non local 
writhings will be inappropriately evaluated by the writhing expression ^ proposed 
by Bouchiat and Mezard, though they make a comparison with the W as measured 
by ([2]). Rossetto and Maggs argue that the Bouchiat and Mezard model works well 
for cases in which the vertical force F is sufficiently high, such that geometrical 
configurations with non local writhing (or winding) lead to a statistically negligible 
contribution. We see in section 4.4| that the Bouchiat and Mezard writhing 
expression actually evaluates a certain proportion of polymer configurations which 
have non-local winding appropriately, that is it returns the same evaluation as the 
polar writhe. 



4-4- A comparison of the various writhing expressions 

We have at several points in this note encountered various writhing expressions. We 
have the double integral writhing expression defined by equation (|2j) which we shall 
denote W, the directional writhing expression used by Bouchiat and Mezard, here 
denoted W z , in their model, finally, we have the polar writhe expression introduced 
above. We use to example curve studies in order to demonstrate that the properties 
of the polar writhe measure are those one would expect of the possible geometrical 
configurations obtainable by a constrained elastic polymer. The procedure will 
involve tracking the evolution of W, W 2 and W p , as applied to a set of curves 
over a period [0, t], as t is increased. 
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Figure 10. AplotofW(x) and W p [x(t)], wherex = (sin Airt, cos 4art,t), evaluated 
over [0, t] for t G [0, 1]. The W p plot is linear, the W becomes linear after t ~ 1.4. 



4-5. Case 1 - Helical configurations 

Helical shapes are a common configuration assumed by symmetric elastic rod or 
polymer, under the force torque pair. Consider a rectilinear elastic rod which central 
axis is place on z, a scenario similar to the constrained DNA experiments we have 
discussed. It is known that a combination of deformations including a force along 
the axis and a torque perpendicular to z will causes buckling modes. One of the 
simpler configurations assume would be a helical shape, that is to say the curves 
axis is helical (see for example Healey and Mehta [17J). Indeed, using the same 
partition function model as Bouchiat and Mezard, Fain et al [14j demonstrated that 
some of the possible energy configurations of the constrained DNA polymer we have 
described can be helical. So the appropriate evaluation of the writhing of helical 
polymer configurations is a desirable properties of our measure. On might imagine 
that the greater the applied torque the greater the possible number of turns of the 
helix. We would expect the writhing to increase linearly with the number of turns 
for a symmetric helix. 



Figure 10 depicts the result of evaluating W and W p of a helix (x(t) = 
(sin 4irt, cos 4nt,t)) over the period [0,1] for t values t e [0,1]. This helix winds 
about a fixed direction (z) and could for example represent the axis of a supercoiled 
DNA molecule bound between two surfaces. We see a marked difference in the 
writhing interpretation of each measure (note that the results for W p in this case 
would be the same as W z as the helical always has a positive vertical gradient). 
Specifically W p increases linearly with t. W increases slower at first until it starts 
to increase linearly after roughly t = 0.4. The reason for this difference derives from 
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(a) 




(c) 




Figure 11. The curve polcv(t), equation) 37), evaluated over a period t £ [— 1, 1], 
for values of a — 0.2 (a), a = 0.8 (b) and a = 1.16 (c) respectively (note the z 
component has been multiplied by a factor of 10 for the sake of aesthetics). In 
(a) the curve has developed a small sigmoidal kink about its midriff. This kink 
develops such that the curve develops yV pn i contributions about its middle section. 
In (c) the curves actually intersect, forming a double point. For a > 1.1.6 the curve 
will have passed through itself leading to a change in its non-local writhing of +2 



(Figure 12 ) 



their interpretation of non-local writhing. W p (x) in this scenario records no non- 
local writhing, over the full parametrization range. W(x) however records non-local 
windings from all possible view points, thus will record the helix as exhibiting non- 
local windings. So we can see that both W z and W p give the appropriate measure 
in such cases, W however does not. To reiterate this is because the helix is a curve 
whose contortion is directionally specific (about z) in this case. The non local writhe 
W is ill equipped to evaluative such curves as much of the directional specificity of 
its geometry will be lost in the averaging process. 
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We now define a curve (depicted in Figure 11), which begins as a curve winding 
locally along z. The coiled section then drops back over its axis leading to non-local 
winding of the curve about itself. It can be parameterized using the following set 
(we call the set polcv, short for polymer-type curve) 

er/(t) = 4= f e-*' 2 dt', (37) 



7T „ ,, 

2\ 







polcv(t) = (a7.5(er/(-2t 2 ) + l)sin27rt 
,2a(er/(-2t 2 ) + l) 

, 5(t + a5erf(-2t 2 ) + l)(t(t - a)(t + a))), (38) 

where a is a constant which controls the contortion of the polymer. Increasing a 
causes a kink to form about the curves midpoint. Note, its endpoints will remain 
fixed. The evolution of this curve has several specific properties which we would 
require our writhe measure to capture, it has a transition from modes which have 
only local writhing about z to those for which distinct sections of the curve are 
wound about each other. As the torque is causing this winding is in the x-y axis 



we would like our measure to reflect this. Finally in (c) Figure 11 the curve passes 
through itself. This the kind of deformation which would be forbidden in the DNA 
experiments. We would like our measure to reflect this fact in some way. 

Figure 12 details the results for W p [x(t)], W z [x(t)] and W[x(£)], as evaluated 
over the period t G [—1,1], for a G [0, 1.4]. There is an agreement in the writhing 
interpretations of W z and W p up until a value of a = 1.16. This includes the period 
in which non-local windings are formed. This occurs as the non local component 
Wpni is such that the sum effectively cancels out W p i components for sections of the 
curve whose gradient along z is negative (see Berger and Prior jl]. This is intriguing 
as it would offer a possible explanation as to why the models of Bouchiat and Mezard 
(and Moroz and Nelson [20]) performed better than one might expect. We note the 
W expression evaluates these curves differently. Once again this is due to the fact 
that the W measure averages aver all directions. 

After a = 1.16 there is a jump in the values of W and W p of +2, as a result 
of the curve passing through itself. W p and W detect the crossing by registering a 
jump of +2. W z on the other hand recognizes only a small difference in writhing 
between two such states. This will be a key factor in the partition function model 
we posit in section |5j The discontinuous jump in writhing due to the curve being 
passed through itself could not be compensated for by a change in twisting, as the 
twisting changes continuously. Thus this change would lead to a change in linking. 

The example studies in this chapter have demonstrated that the polar 
writhe expression satisfies the requirements of key geometrical aspects of possible 
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Figure 12. Plots of W, W z and W p against t for the curve polcv(t), evaluated 
over a period t 6 [—1, 1], over a range of a values a € [0, 1.4]. For a « [0, 1.175] 
Wp and W z agree in their evaluations. However at a ~ 1.75 there is a jump in 
value of both VV and W p as the curve crosses itself. This jump does not occur for 
W z . 



constrained polymer configurations, especially with regards to the DNA micro- 
manipulation experiments in which we are interested. In particular we have seen it 
appropriately captures the directional specificity of the forces applied to this system 
and the self avoidance condition inherent to the molecules. 



5. Linking and the partition function 

The general form of the partition function for an unconstrained molecule would be. 

Z = J V(9,<i),ip) e^, (39) 

Where T>{9 ) cf) ) ip) is the functional space over all possible paths. However we must 
place certain restrictions on the allowed configurations taken by the molecule. As 
discussed in section section [3] this constraint is the net winding. This can be split 
into the local linking of sections for which i = j and non local sections where i ^ j. 
This constraint is applied using the dirac delta function as in [3] . 

Z Z = J V(9,M)5{Z-Ci-Craj e^, (40) 

The partition function has two inputs the applied winding £ and the applied 
stretching force F. It was shown in Berger and Prior that the local linking 
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Figure 13. Figure (a) represents polcv(t) (37), plotted over a range t 6 [— 1, 1], 
with a = 1.1. It represents an example of the curve before it has passed through 
itself (we would expect W p and W z to return the same measure, and they do 
(Figure 12). (b) is the same curve with a = 1.2; the curve has passed through 
itself. 



contribution can be split into twisting and writhing components with the writhing 
component equal to the local polar writhe W p i and the twist given as the integral 
over the angular rotation lo\ ^ (divided by 2ir). For a ribbon i?(x, y), parametrized 
by z, with n turning points about x we have 



T(R, z min , z max ) = — y~] / Oi ( — + I cos 9\ -J ) dz. (11) 



So 



^rriini ^max) T\Ri ^mini ^max) Wplip^j ^mini ^max) ('^) 

We now make a final observation. Consider the energy expression which is an integral 
over the densities defined in section 12.41 

E = E / 1 E 'bend^) + EliM + Ktretch(z) ^ (43) 

J y mi n 

i=l i 

All components of this expression are local. The non local winding only contributes 
to the energy by limiting the set of allowed configurations the polymer can exhibit. 
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Thus we can separate the integral, performing a local partition function for a fixed 
amount of non-local winding (£„/), 



Z{C,F,C 



nil 



^)5((£-C nl )-W pl -r)e^ 



(44) 



One can then integrate over the range of C n i values (note this is theoretically 
unbound) to give us the following. 



Z(C,F) 



Z (£, F, £ n i)d£ n i. 



(45) 



The pair of equations (44) and (45) define a self avoiding partition function. 



6. Conclusions 



We have discussed the geometrical modelling of physical structures which can be 
represented mathematically as ribbon's or rod's, and which are constrained at their 
endpoints. Further they are subject to the condition that all deformations are 
restricted such that they are contained between two planes at the curves endpoints. 
In particular the following results were detailed. 

• We defined the correct components of the ribbon's energy expression, in terms 
of an orthonormal framing of the axis of the ribbon. This differs form previous 
work ( [H] , [5] , [20] ) in that it correctly defined the rates for sections moving 
either upwards (along z) and downwards. The downwards moving sections had 
previously been evaluated incorrectly. Despite this the energy expression was 
shown to be the same. See section [2741 

• We demonstrated that the net winding was the correct topological constraint 
required to restrict the allowed configurations of the ribbon. This expression 
can be defined as as decomposition of local and non-local components. See 
section [3] 

• We derived the correct directional writhing expression for constrained ribbons. 
This derivation used the correct form of a transported magnetic monopole 
(or unit of magnetic flux), this expression is termed the polar writhe (it has 
been derived before in a different fashion |4J). This derivation neatly links 
the problem of a ribbon, manipulated by directionally inclined torsional and 
stretching forces, to the quantum problem of a top in an external magnetic 
field. This builds upon the insight of Bouchiat and Mezard, who first suggested 
this link, by detailing the correct formulation of the monopole filed. Using 
example curve studies we demonstrated that the polar writhe has the required 
geometrical properties for the modelling of constrained ribbons. See section |4} 
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• Using the above results we derive a partition function model for the constrained 
ribbon which forces the constraint of self avoidance using the net winding 
measure. Further, using the decomposition of the net winding into local 
and non-local components, we modified this model such that it allows for the 
possibility of an analytic treatment, see section [5j 

This work represents an improvement to the previous work, discussed in the 
introduction, in that it includes the necessary physical constraint of self- avoidance. 
The next step would be to approach this model using analytical and numerical 



techniques. This could involve minimization of the energy expression (43), subject 
to the net winding constraint, in a similar manner to that detailed by Fain et al 
[14"] . Once could also follow the leads of Bouchiat and Mezard j5] and Moroz and 
Nelson |20j in tackling the partition function. 

Appendix A. Problems with the symmetric top model and twisting 

A key feature of the work of Bouchiat and Mezard is the mapping of the polymer 
chain problem to the quantum mechanical model of a symmetric top, subjected to 
a magnetic field, through an imaginary time transformation (effectively a reversal 
of the Wick rotation). If we assume the molecule is symmetric about its axis then 
the molecule can be modelled as a ribbon which consists of a curve x, representing 
the axis of the molecule and a surrounding curve y which can be used to represent 



one of the phosphate back bones of the DNA molecule (see section Al). In order to 
define the geometry of this molecule we must consider both the contortion (bending 
and winding) of x, in other words the writhing, and the extent to which y wraps 
around x, the twisting. 

We can link this structure and the a top using Euler angle rotations to see 
why the quantum analogy used by Bouchiat and Mezard had some success. Let us 



consider a symmetric top (see (a) in figure A2), consider it shrunk to an infinitesimal 



size with a vector v attached to the molecule and lying in the x — y plane, this vector 
is fixed and will rotate and with the top. This construction is placed at the lowest 



point of the molecule (figure Al ) such that the tip of u lies on y. We now move 
the top along the molecule applying the necessary Euler rotations such that its tip 
points along the tangent direction of x and u will be rotated such that it remains 



on y (figure Al ). Using this set up we can define the geometry of the curve in terms 
of the unique set of Euler angle rotations applied to the top. There is, however, a 
problem with this set-up. 

The issue relates to the twisting of the molecule and the top. The twisting of 
the molecule can be quantified by tracking the rotations, in the x-y plane, of the 
vector u. The key issue is with the use of the Euler rotation 9. The symmetric top 
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Figure Al. Depicted is the process by which an infinitesimal top with an attached 
vector can be used to draw out a polymer consisting of an axis x and a surrounding 
curve y which twists about the axis. This set up could represent a constrained 
DNA molecule. The top begins oriented along the tangent direction of x in the 
plane z — z m im h is then moved up the axis with the necessary Euler rotations 
applied in order for it to be aligned along the tangent direction of x. In addition 
Euler rotations are applied such that the joining vector u stays attached to y. 
This set of Euler rotations can be used to evaluate the twisting of y about x. 



is by definition symmetric about its centre line, which in our set up always lies along 
the tangent direction of x. Imagine spinning the top about its axis when its tip is 
pointing directly upwards and spinning it about its axis in a right handed manner 
(i.e. it follows the right hand thumb rule) (see (a) in figure A2). We ask an observer 
to note the orientation with which it spins. We then (unknown to the observer) 
turn the top such that the tip which was pointing upwards now points downwards 
(in Euler rotations this would require a rotation 9 — ir). To the observer it would 
appear that the top is spinning in the opposite direction (see (b) figure A2). In the 
case of the top representing the polymer molecule this would be a problem. No new 
rotation about z has been applied so the twisting (rotation of the top about its axis 
of symmetry) of the molecule should be the same. However in the symmetric top 
case the twisting appears to have reversed. We need a top which is in some way 
marked (as in (c) figure A2) such, such that we can tell the top from the bottom. 




Figure A2. A depiction of the failure of the symmetric top to distinguish the 
appropriate handedness of twisting. Figure (a) is a symmetric top it has had 
a rotation applied to it in the x-y plane such that it is spinning in the direction 
shown. The top depicted in (b) represents the top in (a) rotated through tt radians 
about the x axis. To the observer it would appear that (b) has the opposite 
rotation to that of (a) . Figure (c) is a top with the same rotation applied to it as 
(a), as shown by the arrows depicting its rotation. However, (c) is not symmetric 
in that its top point has a marker attached to it, making it distinct. We rotate (c) 
through 7r radians (the same rotation as (a) to (b)) in order to produce (d). We 
can now identify that (c) and (d) have the same sense of rotation (right handed 
about the direction point by our marker). 



In this case the observer would be able to note that the top described above has the 
same spin direction after the polar rotation has been applied (compare (c) to (d) in 



figure A2). In terms of Euler angle rotation this is achieved by using | cos#| rather 
than cos 9 in our calculations. The modulus sign will mean that there is a symmetry 
in the twisting of the top about 9 = 71/ 2. 
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